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ABSTRACT 

We present a pedagogical and phenomenological introduction to the study of cosmic 
microwave background (CMB) polarization to build intuition about the prospects and 
challenges facing its detection. Thomson scattering of temperature anisotropies on the 
last scattering surface generates a linear polarization pattern on the sky that can be 
simply read off from their quadrupole moments. These in turn correspond directly 
to the fundamental scalar (compressional), vector (vortical), and tensor (gravitational 
wave) modes of cosmological perturbations. We explain the origin and phenomenology 
of the geometric distinction between these patterns in terms of the so-called electric 
and magnetic parity modes, as well as their correlation with the temperature pattern. 
By its isolation of the last scattering surface and the various perturbation modes, the 
polarization provides unique information for the phenomenological reconstruction of the 
cosmological model. Finally we comment on the comparison of theory with experimental 
data and prospects for the future detection of CMB polarization. 
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1. Introduction 

Why should we be concerned with the polarization 
of the cosmic microwave background (CMB) anisotro¬ 
pies? That the CMB anisotropies are polarized is a 
fundamental prediction of the gravitational instability 
paradigm. Under this paradigm, small fluctuations in 
the early universe grow into the large scale structure we 
see today. If the temperature anisotropies we observe 
are indeed the result of primordial fluctuations, their 
presence at last scattering would polarize the CMB an¬ 
isotropies themselves. The verification of the (partial) 
polarization of the CMB on small scales would thus rep¬ 
resent a fundamental check on our basic assumptions 
about the behavior of fluctuations in the universe, in 
much the same way that the redshift dependence of the 
CMB temperature is a test of our assumptions about 
the background cosmology. 

Furthermore, observations of polarization provide an 
important tool for reconstructing the model of the fluc¬ 
tuations from the observed power spectrum (as distinct 
from fitting an a priori model prediction to the observa¬ 
tions) . The polarization probes the epoch of last scat¬ 
tering directly as opposed to the temperature fluctua¬ 
tions which may evolve between last scattering and the 
present. This localization in time is a very powerful 
constraint for reconstructing the sources of anisotropy. 
Moreover, different sources of temperature anisotropies 
(scalar, vector and tensor) give different patterns in 
the polarization: both in its intrinsic structure and in 
its correlation with the temperature fluctuations them¬ 
selves. Thus by including polarization information, one 
can distinguish the ingredients which go to make up the 
temperature power spectrum and so the cosmological 
model. 

Finally, the polarization power spectrum provides 
information complementary to the temperature power 
spectrum even for ordinary (scalar or density) pertur¬ 
bations. This can be of use in breaking parameter de¬ 
generacies and thus constraining cosmological parame¬ 
ters more accurately. The prime example of this is the 
degeneracy, within the limitations of cosmic variance, 
between a change in the normalization and an epoch of 
“late” reionization. 

Yet how polarized are the fluctuations? The degree 
of linear polarization is directly related to the quadru- 
pole anisotropy in the photons when they last scatter. 
While the exact properties of the polarization depend 
on the mechanism for producing the anisotropy, several 
general properties arise. The polarization peaks at an¬ 


gular scales smaller than the horizon at last scattering 
due to causality. Furthermore, the polarized fraction 
of the temperature anisotropy is small since only those 
photons that last scattered in an optically thin region 
could have possessed a quadrupole anisotropy. The frac¬ 
tion depends on the duration of last scattering. For the 
standard thermal history, it is 10% on a characteristic 
scale of tens of arcminutes. Since temperature aniso¬ 
tropies are at the 10~ 5 level, the polarized signal is at 
(or below) the 10~ 6 level, or several /rK, representing a 
significant experimental challenge. 


Our goal here is to provide physical intuition for these 
issues. For mathematical details, we refer the reader 


to Kamionkowski et al. (1997), Zaldarriaga & Seljak 


(1 997), |H u fe W hite| ( 1997) as w ell as pion eering work 
by Bond fc Efstathiou (1984) and Polnarev (1985). The 
outline of the paper is as follows. We begin in §|] by 
examining the general properties of polarization forma¬ 
tion from Thomson scattering of scalar, vector and ten¬ 
sor anisotropy sources. We discuss the properties of 
the resultant polarization patterns on the sky and their 
correlation with temperature patterns in (j ; 3|. These con¬ 
siderations are applied to the reconstruction problem in 
§[f|. The current state of observations and techniques for 
data analysis are reviewed in §||. We conclude in §||with 
comments on future prospects for the measurement of 
CMB polarization. 


2. Thomson Scattering 


2.1. From Anisotropies to Polarization 


The Thomson scattering cross section depends on po¬ 


larization as (see e.g. Chandrasekhar 1960) 


dar 

dfl 


( 1 ) 


where e (e') are the incident (scattered) polarization di¬ 
rections. Heuristically, the incident light sets up oscilla¬ 
tions of the target electron in the direction of the elec¬ 
tric field vector E, i.e. the polarization. The scattered 
radiation intensity thus peaks in the direction normal 
to, with polarization parallel to, the incident polariza¬ 
tion. More formally, the polarization dependence of the 
cross section is dictated by electromagnetic gauge in¬ 
variance and thus follows from very basic principles of 
fundamental physics. 

If the incoming radiation field were isotropic, orthog¬ 
onal polarization states from incident directions sepa¬ 
rated by 90° would balance so that the outgoing radia¬ 
tion would remain unpolarized. Conversely, if the inci- 
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Fig. 1.— Thomson scattering of radiation with a quad- 
rupole anisotropy generates linear polarization. Blue 
colors (thick lines) represent hot and red colors (thin 
lines) cold radiation. 

dent radiation field possesses a quadrupolar variation in 
intensity or temperature (which possess intensity peaks 
at 90° = 7t/2 separations), the result is a linear polar¬ 
ization of the scattered radiation (see Fig. jl|). A reversal 
in sign of the temperature fluctuation corresponds to a 
90° rotation of the polarization, which reflects the spin- 
2 nature of polarization. 

In terms of a multipole decomposition of the ra¬ 
diation field into spherical harmonics, Y f m (0, (f>), the 
five quadrupole moments are represented by l = 2, 
m = 0, ±1, ±2. The orthogonality of the spherical har¬ 
monics guarantees that no other moment can generate 
polarization from Thomson scattering. In these spheri¬ 
cal coordinates, with the north pole at 9 = 0, we call a 
N-S (E-W) polarization component Q > 0 (Q < 0) and 
a NE-SW (NW-SE) component U > 0 (U < 0). The 
polarization amplitude and angle clockwise from north 
are 

P = VQ 2 + u 2 , a=itan ~ 1 {U/Q). (2) 

Alternatively, the Stokes parameters Q and U repre¬ 
sent the diagonal and off diagonal components of the 
symmetric, traceless, 2x2 intensity matrix in the po¬ 


larization plane spanned by (eg, e^), 

E*Ej - i. SijE 2 a Qa 3 + Ua 1 , (3) 

where cq are the Pauli matrices and circular polarization 
is assumed absent. 

If Thomson scattering is rapid, then the randomiza¬ 
tion of photon directions that results destroys any quad¬ 
rupole anisotropy and polarization. The problem of un¬ 
derstanding the polarization pattern of the CMB thus 
reduces to understanding the quadrupolar temperature 
fluctuations at last scattering. 

Temperature perturbations have 3 geometrically dis¬ 
tinct sources: the scalar (compressional), vector (vorti¬ 
cal) and tensor (gravitational wave) perturbations. For¬ 
mally, they form the irreducible basis of the symmetric 
metric tensor. We shall consider each of these below 
and show that the scalar, vector, and tensor quadru¬ 
pole anisotropy correspond to m = 0,±1,±2 respec¬ 
tively. This leads to different patterns of polarization 
for the three sources as we shall discuss in §j|. 



Fig. 2.— The scalar quadrupole moment {£ = 2,m = 
0). Flows from hot (blue) regions into cold (red), v || k, 
produce the azimuthally symmetric pattern Y 2 ° depicted 
here. 

2.2. Scalar Perturbations 

The most commonly considered and familiar types of 
perturbations are scalar modes. These modes represent 


3 












Fig. 3.— The transformation of quadrupole anisotropies into linear polarization, (a) The orientation of the quadru- 
pole moment with respect to the scattering direction n determines the sense and magnitude of the polarization. It 
is aligned with the cold (red, long) lobe in the eg®e,p tangent plane, (b) In spherical coordinates where nk = cos 9, 
the polarization points north-south (Q) with magnitude varying as sin 2 9 for scalar fluctuations. 


perturbations in the (energy) density of the cosmolog¬ 
ical fluid(s) at last scattering and are the only fluctu¬ 
ations which can form structure though gravitational 
instability. 

Consider a single large-scale Fourier component of 
the fluctuation, i.e. for the photons, a single plane wave 
in the temperature perturbation. Over time, the tem¬ 
perature and gravitational potential gradients cause a 
bulk flow, or dipole anisotropy, of the photons. Both 
effects can be described by introducing an “effective” 
temperature 


(AT/T) eff = A T/T + T , (4) 

where is the gravitational potential. Gradients in 
the effective temperature always create flows from hot 
to cold effective temperature. Formally, both pressure 
and gravity act as sources of the momentum density of 
the fluid in a combination that is exactly the effective 
temperature for a relativistic fluid. 

To avoid confusion, let us explicitly consider the case 
of adiabatic fluctuations, where initial perturbations to 
the density imply potential fluctuations that dominate 
at large scales. Here gravity overwhelms pressure in 
overdense regions causing matter to flow towards den¬ 
sity peaks initially. Nonetheless, overdense regions are 
effectively cold initially because photons must climb out 
of the potential wells they create and hence lose energy 
in the process. Though flows are established from cold 
to hot temperature regions on large scales, they still 


go from hot to cold effective temperature regions. This 
property is true more generally of our adiabatic assump¬ 
tion: we hereafter refer only to effective temperatures 
to keep the argument general. 

Let us consider the quadrupole component of the 
temperature pattern seen by an observer located in a 
trough of a plane wave. The azimuthal symmetry in 
the problem requires that v || k and hence the flow is 
irrotational V x v = 0. Because hotter photons from 
the crests flow into the trough from the ±k directions 
while cold photons surround the observer in the plane, 
the quadrupole pattern seen in a trough has an m = 0, 

Y° oc 3 cos 2 9 — 1, (5) 


structure with angle n k = cos 9 (see Fig. |2). The oppo¬ 
site effect occurs at the crests, reversing the sign of the 
quadrupole but preserving the m = 0 nature in its local 
angular dependence. The full effect is thus described 
by a local quadrupole modulated by a plane wave in 
space, — y 2 °(n) exp(?'k • x), where the sign denotes the 
fact that photons flowing into cold regions are hot. This 
infall picture must be modified slightly on scales smaller 
than the sound horizon where pressure plays a role (see 
§3.3.2), however the essential property that the flows are 
parallel to k and thus generate an m = 0 quadrupole 
remains true. 


The sense of the quadrupole moment determines the 
polarization pattern through Thomson scattering. Re¬ 
call that polarized scattering peaks when the tempera- 
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Fig. 4.- Polarization pattern for £ = 2, to = 0, note the azimuthal symmetry. The scattering of a scalar 
to = 0 quadrupole perturbation generates the electric E (yellow, thick lines) pattern on the sphere. Its ro¬ 
tation by 45° represents the orthogonal magnetic B (purple, thin lines) pattern. Animation (available at 
ittp://www. sns. ias .edu/~whu/polar/scalaran.html): as the line of sight ri changes, the lobes of the quadru¬ 
pole rotate in and out of the tangent plane. The polarization follows the orientation of the colder (red) lobe in the 
tangent plane. 


ture varies in the direction orthogonal to n. Consider 
then the tangent plane eg (g) with normal n. This 
may be visualized in an angular “lobe” diagram such 
as Fig. U as a plane which passes through the “origin” 
of the quadrupole pattern perpendicular to the line of 
sight. The polarization is maximal when the hot and 
cold lobes of the quadrupole are in this tangent plane, 
and is aligned with the component of the colder lobe 
which lies in the plane. As 9 varies from 0 to n/2 (pole 
to equator) the temperature differences in this plane in¬ 
crease from zero (see Fig. So- The local polarization at 
the crest of the temperature perturbation is thus purely 
in the N-S direction tapering off in amplitude toward 
the poles (see Fig. gb). This pattern represents a pure 
Q-field on the sky whose amplitude varies in angle as 
an £ = 2, m = 0 tensor or spin-2 spherical harmonic 

Q = sin 2 9, U = 0. (6) 

In different regions of space, the plane wave modulation 
of the quadrupole can change the sign of the polariza¬ 
tion, but not its sense. 

This pattern (Fig. g, yellow lines) is of course not the 
only logical possibility for an £ = 2, to = 0 polarization 
pattern. Its rotation by 45° is also a valid configuration 
(purple lines). This represents a pure NW-SE (and by 
sign reversal NE-SW), or [/-polarization pattern. We 


return in §3.3 to consider the geometrical distinction 
between the two patterns, the electric and magnetic 
modes. 



Fig. 5.— The vector quadrupole moment (£ = 2, 
to = 1). Since v _L k, the Doppler effect generates a 
quadrupole pattern with lobes 45° from v and k that is 
spatially out of phase (interplane peaks) with v. 
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Fig. 6. - Polarization pattern for £ = 2, m = 1. The scattering of a vector (m = 1) quadrupole perturbation 
generates th e E pattern (yellow, thick lines) as opposed to the B, (purple, thin lines) pattern. Animation 
(available at ittp://www.sns.ias.edu/^whu/polar/vectoran.html): same as for scalars. 


2.3. Vector Perturbations 

Vector perturbations represent vortical motions of 
the matter, where the velocity field v obeys V ■ v = 0 
and V x v / 0, similar to “eddies” in water. There is 
no associated density perturbation, and the vorticity is 
damped by the expansion of the universe as are all mo¬ 
tions that are not enhanced by gravity. However, the 
associated temperature fluctuations, once generated, do 
not decay as both AT and T scale similarly with the 
expansion. For a plane wave perturbation, the veloc¬ 
ity field v _L k with direction reversing in crests and 
troughs (see Fig. ||). The radiation field at these ex¬ 
trema possesses a dipole pattern due to the Doppler 
shift from the bulk motion. Quadrupole variations van¬ 
ish here but peak between velocity extrema. To see this, 
imagine sitting between crests and troughs. Looking up 
toward the trough, one sees the dipole pattern projected 
as a hot and cold spot across the zenith; looking down 
toward the crest, one sees the projected dipole reversed. 
The net effect is a quadrupole pattern in temperature 
with m = ±1 

Y^ 1 ocsindcos9e ±l ^. (7) 

The lobes are oriented at 45° from k and v since the line 
of sight velocity vanishes along k and at 90 degrees to k 
here. The latter follows since midway between the crests 
and troughs v itself is zero. The full quadrupole distri¬ 
bution is therefore described by —iV 2 ±:L (n) exp(zk • x), 
where i represents the spatial phase shift of the quad¬ 


rupole with respect to the velocity. 

Thomson scattering transforms the quadrupole tem¬ 
perature anisotropy into a local polarization field as be¬ 
fore. Again, the pattern may be visualized from the 
intersection of the tangent plane to n with the lobe pat¬ 
tern of Fig. ||. At the equator (0 = 7r/2), the lobes are 
oriented 45° from the line of sight n and rotate into and 
out of the tangent plane with </>. The polarization pat¬ 
tern here is a pure [/-field which varies in magnitude as 
sin <f>. At the pole 0 = 0, there are no temperature vari¬ 
ations in the tangent plane so the polarization vanishes. 
Other angles can equally well be visualized by viewing 
the quadrupole pattern at different orientations given 
by n. 

The full £ = 2, m = 1 pattern, 

Q = -sin0cos0e®^, U = — ismOe 1 ^ (8) 

is displayed explicitly in Fig. || (yellow lines, real part). 
Note that the pattern is dominated by [/-contributions 
especially near the equator. The similarities and dif¬ 
ferences with the scalar pattern will be discussed more 
fully in §|| 

2.4. Tensor Perturbations 

Tensor fluctuations are transverse-traceless pertur¬ 
bations to the metric, which can be viewed as gravita¬ 
tional waves. A plane gravitational wave perturbation 
represents a quadrupolar “stretching” of space in the 
plane of the perturbation (see Fig. fj]). As the wave 
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3. Polarization Patterns 



Tensors 
(Gravity Waves) 


Fig. 7. The tensor quadrupole moment (m = 2). 
Since gravity waves distort space in the plane of the 
perturbation, changing a circle of test particles into an 
ellipse, the radiation acquires an to = 2 quadrupole 
moment. 


passes or its amplitude changes, a circle of test par¬ 
ticles in the plane is distorted into an ellipse whose 
semi-major axis —> semi-minor axis as the spatial phase 
changes from crest —► trough (see Fig. |l], yellow ellipses). 
Heuristically, the accompanying stretching of the wave¬ 
length of photons produces a quadrupolar temperature 
variation with an to = ±2 pattern 

Y 2 ±2 oc sin 2 0e ±2i?i (9) 

in the coordinates defined by k. 

Thomson scattering again produces a polarization 
pattern from the quadrupole anisotropy. At the equa¬ 
tor, the quadrupole pattern intersects the tangent (e^0 
e^) plane with hot and cold lobes rotating in and out 
of the direction with the azimuthal angle <j >. The po¬ 
larization pattern is therefore purely Q with a cos(2(f>) 
dependence. At the pole, the quadrupole lobes lie com¬ 
pletely in the polarization plane and produces the maxi¬ 
mal polarization unlike the scalar and vector cases. The 
full pattern, 

Q = (1 + cos 2 0)e 2l( ^, U = — 2icos9e 21 ^, (10) 


The considerations of (ji| imply that scalars, vec¬ 
tors, and tensors generate distinct patterns in the polar¬ 
ization of the CMB. However, although they separate 
cleanly into to = 0, ±1,±2 polarization patterns for a 
single plane wave perturbation in the coordinate sys¬ 
tem referenced to k, in general there will exist a spec¬ 
trum of fluctuations each with a different k. Therefore 
the polarization pattern on the sky does not separate 
into m = 0, ±1, ±2 modes. In fact, assuming statistical 
isotropy, one expects the ensemble averaged power for 
each multipole t to be independent of to. Nonetheless, 
certain properties of the polarization patterns discussed 
in the last section do survive superposition of the per¬ 
turbations: in particular, its parity and its correlation 
with the temperature fluctuations. We now discuss how 
one can describe polarization patterns on the sky arising 
from a spectrum of k modes. 

3.1. Electric and Magnetic Modes 

Any polarization pattern on the sky can be separated 
into “electric” (E) and “magnetic” (E) components^. 
This decomposition is useful both observationally and 
theoretically, as we will discuss below. There are two 
equivalent ways of viewing the modes that reflect their 
global and local properties respectively. The nomencla¬ 
ture reflects the global property. Like multipole radi¬ 
ation, the harmonics of an E-mode have (— l) e parity 
on the sphere, whereas those of a E-mode have (—1 Y +1 
parity. Under n —> —n, the E-mode thus remains un¬ 
changed for even £, whereas the E-mode changes sign as 
illustrated for the simplest case l = 2, to = 0 in Fig. || 
(recall that a rotation by 90° represents a change in 
sign). Note that the E and E multipole patterns are 
45° rotations of each other, i.e. Q —> U and U —> —Q. 
Since this parity property is obviously rotationally in¬ 
variant, it will survive integration over k. 

The local view of E and E-modes involves the sec¬ 
ond derivatives of the polarization amplitude ( second 
derivatives because polarization is a tensor or spin-2 
object). In much the same way that the distinction be¬ 
tween electric and magnetic fields in electromagnetism 
involves vanishing of gradients or curls (i.e. first deriva¬ 
tives) for the polarization there are conditions on the 
second (covariant) derivatives of Q and U. For an E- 
mode, the difference in second (covariant) derivatives 


is shown in Fig. g (real part). Note that Q and U are 
present in nearly equal amounts for the tensors. 


2 These com ponents are called the 


ponents by Kamionkowski et al 


“grad” (G) and “curl” (C) com- 
( 1997 ). 
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Fig. 8.- Polarization pattern for l = 2 m = 2. Scattering of a tensor to = 2 perturbation generates the 
E (yellow, thick lines) pattern as opposed to the E (purple, thin lines) pattern. Animation (available at 
|ittp: //www. sns. ias . edu/^whu/polar/tensoran.html): same as for scalars and vectors. 



Fig. 9.— The electric (E) and magnetic ( B ) modes are distinguished by their behavior under a parity transfor¬ 
mation n —> —n. E-modes have (—1)^ and E-modes have (—1)^ +1 parity; here (£ = 2, to = 0), even and odd 
respectively. The local distinction between the two is that the polarization axis is aligned with the principle axes 
of the polarization amplitude for E and crossed with them for B. Dotted lines represent a sign reversal in the 
polarization. 


of U along eg and vanishes as does that for Q along 
eg + and eg — e^. For a E-mode, Q and U are in¬ 
terchanged. Recalling that a Q-field points in the eg or 
direction and a E-field in the crossed direction, we 
see that the Hessian or curvature matrix of the polar¬ 
ization amplitude has principle axes in the same sense 


as the polarization for E and 45° crossed with it for E 
(see Fig. |j). Stated another way, near a maximum of 
the polarization (where the first derivative vanishes) the 
direction of greatest change in the polarization is paral¬ 
lel/perpendicular and at 45° degrees to the polarization 
in the two cases. 











The distinction is best illustrated with examples. 
Take the simplest case of £ = 2, to = 0 where the 
E-mode is a Q = sin 2 9 field and the B-mode is a 
U = sin 2 9 field (see Fig. |^). In both cases, the ma¬ 
jor axis of the curvature lies in the eg direction. For 
the E-mode, this is in the same sense; for the B-mode 
it is crossed with the polarization direction. The same 
holds true for the m = 1,2 modes as can be seen by 
inspection of Fig. and || 


3.2. Electric and Magnetic Spectra 


Thomson scattering can only produce an E-mode lo¬ 
cally since the spherical harmonics that describe the 
temperature anisotropy have (—1) £ electric parity. In 
Figs. [|, ||, and the £ = 2, m = 0,1, 2 E-mode patterns 
are shown in yellow. The B-mode represents these pat¬ 
terns rotated by 45° and are shown in purple and cannot 
be generated by scattering. In this way, the scalars, vec¬ 
tors, and tensors are similar in that scattering produces 
a local £ = 2 E-mode only. 

However, the pattern of polarization on the sky is 
not simply this local signature from scattering but is 
modulated over the last scattering surface by the plane 
wave spatial dependence of the perturbation (compare 
Figs. |] and [Io|) . The modulation changes the amplitude, 
sign, and angular structure of the polarization but not 
its nature, e.g. a Q-polarization remains Q. Nonethe¬ 
less, this modulation generates a B-mode from the local 
E-mode pattern. 

The reason why this occurs is best seen from the 
local distinction between E and B-modes. Recall that 
E-modes have polarization amplitudes that change par¬ 
allel or perpendicular to, and B-modes in directions 45° 
away from, the polarization direction. On the other 
hand, plane wave modulation always changes the po¬ 
larization amplitude in the direction k or N-S on the 
sphere. Whether the resultant pattern possesses E or 
B-contributions depends on whether the local polariza¬ 
tion has Q or [/-contributions. 


For scalars, the modulation is of a pure Q-field and 
thus its E-mode nature is preserved (Kamionkowski et 
Fl| 1997; Zaldarriaga & Seljak 1997). For the vectors, 
the U -mode dominates the pattern and the modulation 
is crossed with the polarization direction. Thus vectors 
generate mainly B-modes for short wavelength fluctua¬ 
tions ( Hu fe White| 1997). For the tensors, the compa¬ 
rable Q and U components of the local pattern imply 
a more comparable distribution of E and B modes at 
short wavelengths (see Fig. 0a). 



Fig. 10.- Modulation of the local pattern Fig. |b by 
plane wave fluctuations on the last scattering surface. 
Yellow points represent polarization out of the plane 
with magnitude proportional to sign. The plane wave 
modulation changes the amplitude and sign of the po¬ 
larization but does not mix Q and U. Modulation can 
mix E and B however if U is also present. 


These qualitative considerations can be quantified by 
noting that plane wave modulation simply represent the 
addition of angular momentum from the plane wave 
(Y/ 3 ) with the local spin angular dependence. The re¬ 
sult is that plane wave modulation takes the £ = 2 lo¬ 
cal angular dependence to higher £ (smaller angles) and 
splits the signal into E and B components with ratios 
which are related to Clebsch-Gordan coefficients. At 
short wavelengths, these ratios are B/E = 0,6,8/13 in 
power for sca lars, vectors, and tensors (see Fig. 0b and 
Hu fe Whitel 1997). 

The distribution of power in multipole [-space is also 
important. Due to projection, a single plane wave con¬ 
tributes to a range of angular scales £ < kr where r 
is the comoving distance to the last scattering surface. 
From Fig. |Io], we see that the smallest angular, largest 
£ ss kr variations occur on lines of sight n • k = 0 or 
9 = 7t/2 though a small amount of power projects to 
£ <C kr as 9 —> 0. The distribution of power in multipole 
space of Fig. 0b can be read directly off the local polar¬ 
ization pattern. In particular, the region near 9 = n/2 
shown in Fig. |Tl| a determines the behavior of the main 
contribution to the polarization power spectrum. 

The full power spectrum is of course obtained by 
summing these plane wave contributions with weights 
dependent on the source of the perturbations and the 
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Fig. 11. - The E and B components of a plane wave perturbation, (a) Modulation of the local E-quadrupole 
pattern (yellow) from scattering by a plane wave. Modulation in the direction of (or orthogonal to) the polarization 
generates an FI-mode with higher l\ modulation in the crossed (45°) direction generates a E-mode with higher i. 
Scalars generate only .E-modes, vectors mainly E-modes, and tensors comparable amounts of both, (b) Distribution 
of power in a single plane wave with kr = 100 in multipole t from the addition of spin and orbital angular momentum. 
Features in the power spectrum can be read directly off the pattern in (a). 


dynamics of their evolution up to last scattering. Sharp 
features in the fc-power spectrum will be preserved in 
the multipole power spectrum to the extent that the 
projectors in Fig. lib approximate delta functions. For 
scalar E-modes, the sharpness of the projection is en¬ 
hanced due to strong Q-contributions near 9 = 7r/2 
{£ ~ kr) that then diminish as 9 —■> 0 (E <C kr). The 
same enhancement occurs to a lesser extent for vector 
E-modes due to U near 7t/2 and tensor E-modes due 
to Q there. On the other hand, a supression occurs for 
vector E and tensor E-modes due to the absence of Q 
and U at 7t/2 respectively. These considerations have 
consequences for the sharpness of features in the polar¬ 
ization power spectrum, and the generation of asymp¬ 
totic “ta ils” to the po larization spectrum at low-f (see 
§4.4 and Hu fe White 1997) . 


3.3. Temperature-Polarization Correlation 

As we have seen in §^, the polarization pattern re¬ 
flects the local quadrupole anisotropy at last scattering. 


Hence the temperature and polarization anisotropy pat¬ 
terns are correlated in a way that can distinguish be¬ 
tween the scalar, vector and tensor sources. 

There are two subtleties involved in establishing the 
correlation. First, the quadrupole moment of the tem¬ 
perature anisotropy at last scattering is not generally 
the dominant source of anisotropies on the sky, so the 
correlation is neither 100% nor necessarily directly vis¬ 
ible as patterns in the map. 

The second subtlety is that the correlation occurs 
through the E-mode unless the polarization has been 
Faraday or otherwise rotated between the last scatter¬ 
ing surface and the present. As we have seen an E- 
mode is modulated in the direction of, or perpendicular 
to, its polarization axis. To be correlated with the tem¬ 
perature, this modulation must also correspond to the 
modulation of the temperature perturbation. The two 
options are that E is parallel or perpendicular to crests 
in the temperature perturbation. As modes of different 
direction k are superimposed, this translates into a ra¬ 
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dial or tangential polarization pattern around hot spots 
(see Fig. |l2|a). 

On the other hand E-modes do not correlate with the 
temperature. In other words, the rotation of the pat¬ 
tern in Fig. [Ti| i by 45° into those of Fig. E§) (solid and 
dashed lines) cannot be generated by Thomson scatter¬ 
ing. The temperature held that generates the polariza¬ 
tion has no way to distinguish between points reflected 
across the symmetric hot spot and so has no way to 
choose between the ±45° rotations. This does not how¬ 
ever imply that B vanishes. For example, for a single 
plane wave fluctuation B can change signs across a hot 
spot and hence preserve reflection symmetry (e.g. Fig. |] 
around the hot spot 9 = 7r/4, <j> = 0). However superpo¬ 
sition of oppositely directed waves as in Fig. |i~2| b would 
destroy the correlation with the hot spot. 

7"hc problem of understanding the correlation thus 


the temperature perturbation on the last scattering sur¬ 
face itself. The Doppler contributions can be up to half 
of the total contribution but as we have seen in §2.3 
do not correlate with the quadrupole moment. Contri¬ 
butions after last scattering, while potentially strong in 
isocurvature models for example, also rapidly lose their 
correlation with the quadrupole at last scattering. 

As we have seen, the temperature gradient associ¬ 
ated with the scalar fluctuation makes the photon fluid 
flow from hot regions to cold initially. Around a point 
on a crest therefore the intensity peaks in the directions 
along the crest and falls off to the neighboring troughs. 
This corresponds to a polarization perpendicular to the 
crest (see Fig. jli|). Around a point on a trough the 
polarization is parallel to the trough. As we superpose 
waves with different k we find the pattern is tangential 


around hot s pots and radial around cold spots (Crit 


bre aks down into two steps: (1) determine how the tenden ct al. 1995). It is important to stress that the 


quadrupole moment of the temperature at last scatter¬ 
ing correlates with the dominant source of anisotropies; 
(2) isolate the .E-component (Q-component in k coordi¬ 
nates) and determine whether it represents polarization 
parallel or perpendicular to crests and so radial or tan¬ 
gential to hot spots. 
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Fig. 12.— Temperature-polarization cross correlation. 
E-parity polarization perpendicular (parallel) to crests 
generates a tangential (radial) polarization held around 
hot spots. E-parity polarization does not correlate with 
temperature since the ±45° rotated contributions from 
oppositely directed modes cancel. 


hot and cold spots refer only to the temperature com¬ 
ponent which is correlated with the polarization. The 
correlation increases at scales approaching the horizon 
at last scattering since the quadrupole anisotropy that 
generates polarization is caused by flows. 

For the vectors, no temperature perturbations exist 
on the last scattering surface and again Doppler contri¬ 
butions do not correlate with the quadrupole. Thus the 
main correlations with the temperature will come from 
the quadrupole moment itself. The correlated signal 
is reduced since the strong E-contributions of vectors 
play no role. Hot spots occur in the direction 9 = 7t/4, 
4> = 0 where the hot lobe of the quadrupole is pointed 
at the observer (see Fig. |j). Here the Q (E) component 
lies in the N-S direction perpendicular to the crest (see 
Fig. |(i]). Thus the pattern is tangential to hot spot, like 
scalars ( |1 In fe White 1997). The signature peaks near 
the horizon at last scattering for reasons similar to the 
scalars. 

For the tensors, both the temperature and polariza¬ 
tion perturbations arise from the quadrupole moment, 
which fixes the sense of the main correlation. Hot 
spots and cold spots occur when the quadrupole lobe 
is pointed at the observer, 9 = 7r/2, (j) = 7r/2 and 37t/2. 
The cold lobe and hence the polarization then points 
in the E-W direction. Unlike the scalars and vectors, 


the pattern will be mainly radial to hot spots (Critten¬ 


den et al. 1995). Again the polarization and hence the 


3.3.1. Large Angle Correlation Pattern 

Consider first the large-angle scalar perturbations. 
Here the dominant source of correlated anisotropies is 


cross-correlation peaks near the horizon at last scatter¬ 
ing since gravitational waves are frozen before horizon 
crossing (Polnarev 1985). 
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Fig. 13.— Time evolution of acoustic oscillations. 
The polarization is related to the flows v (red) which 
form quadrupole anisotropies such that its product 
(green) with the effective temperature (blue) reflects 
the temperature-polarization cross correlation. As de¬ 
scribed in the text the adiabatic (a) and isocurvature 
(b) modes differ in the phase of the oscillation in all 
three quantities. Temperature and polarization are an¬ 
ticorrelated in both cases at early times or large scales 
ks < 1. 


3.3.2. Small Angle Correlation Pattern 

Until now we have implicitly assumed that the evolu¬ 
tion of the perturbations plays a small role as is gener¬ 
ally true for scales larger than the horizon at last scat¬ 
tering. Evolution plays an important role for small- 
scale scalar perturbations where there is enough time 
for sound to cross the perturbation before last scatter¬ 
ing. The infall of the photon fluid into troughs com¬ 
presses the fluid, increasing its density and temperature. 
For adiabatic fluctuations, this compression reverses the 
sign of the effective temperature perturbation when the 
sound horizon s grows to be ks ~ n/2 (see Fig. |Tg|a). 
This reverses the sign of the correlation with the quad¬ 
rupole moment. Infall continues until the compression 
is so great that photon pressure reverses the flow when 


ks « 7T. Again the correlation reverses sign. This pat¬ 
tern of correlations and anticorrelations continues at 
twice the frequency of the acoustic oscillations them¬ 
selves (see Fig. &). Of course the polarization is only 
generated at last scattering so the correlations and an¬ 
ticorrelations are a function of scale with sign changes 
at multiples of tt/2 s*, where s* is the sound horizon at 
last scattering. As discussed in ; ; [l.2| . these fluctuations 
project onto anisotropies as £ ~ kr. 


Any scalar fluctuation will obey a similar pattern 
that reflects the acoustic motions of the photon fluid. 
In particular, at the largest scales the fcs* <C 7r/2, the 
polarization must be anticorrelated with the tempera¬ 
ture because the fluid will always flow with the tem¬ 
perature gradient initially from hot to cold. However, 
where the sign reversals occur depend on the acoustic 
dynamics and so is a useful probe of the nature of the 
scalar perturb ations, e.g. w hether they are adiabatic or 
isocurvature ( Hu fe White 1996). In typical isocurva¬ 
ture models, the lack of initial temperature perturba¬ 
tions delays the acoustic oscillation by 7r/2 in phase so 
that correlations reverse at /cs* = 7r (see Fig. &)• 


For the vector and tensor modes, strong evolution 
can introduce a small correlation with temperature fluc¬ 
tuations generated after last scattering. The effect is 
generally weak and model dependent and so we shall 
not consider it further here. 


4. Model Reconstruction 

While it is clear how to compare theoretical predic¬ 
tions of a given model with observations, the recon¬ 
struction of a phenomenological model from the data 
is a more subtle issue. The basic problem is that in 
the CMB, we see the whole history of the evolution in 
redshift projected onto the two dimensional sky. The re¬ 
construction of the evolutionary history of the universe 
might thus seem an ill-posed problem. 

Fortunately, one needs only to assume the very basic 
properties of the cosmological model and the gravita¬ 
tional instability picture before useful information may 
be extracted. The simplest example is the combina¬ 
tion of the amplitude of the temperature fluctuations, 
which reflect the conditions at horizon crossing, and 
large scale structure today. Another example is the 
acoustic peaks in the temperature which form a snap¬ 
shot of conditions at last scattering on scales below the 
horizon at that time. In most models, the acoustic sig¬ 
nature provides a wealth of information 
ical parameters and structure formation 


on cosmolog- 
(Hu fe White 
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Fig. 14.— Temperature, polarization, and temperature-polarization cross correlation predictions and sensitivity 
of MAP for a fiducial model flo = 1, 12s = 0.1, h = 0.5 cold dark matter. The raw MAP satellite sensitivity 
(1 a errors on the recovered power spectrum binned in £) is approximated by noise weights of w^ 1 = (O.ll^tA') 2 
for the temperature and Wp 1 = (0.15/zAT) 2 for the polarization and a FWHM beam of 0.25°. Note that errors 
between the spectra are correlated. While the reionized model (purple r = 0.1) is impossible to distinguish 
from the fiducial model from temperature anisotropies alone, its effect on polarization is clearly visible at low 
l. Dashed lines for the temperature-polarization corr elation represent anticorrelation. Animation (available at 
ittp: //www. sns. ias . edu/^whu/polar/tauan .html): Variations in the spectra as r is stepped from 0—1. 
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1996). Unfortunately, it does not directly tell us the 
behavior on the largest scales where important causal 
distinctions between models lie. Furthermore it may be 
absent in models with complex evolution on small scales 
such as cosmological defect models. 

Here we shall consider how polarization information 
aids the reconstruction process by isolating the last scat¬ 
tering surface on large scales, and separating scalar, vec¬ 
tor and tensor components. If and when these proper¬ 
ties are determined, it will become possible to establish 
observationally the basic properties of the cosmological 
model such as the nature of the initial fluctuations, the 
mechanism of their generation, and the thermal history 
of the universe. 

4.1. Last Scattering 


4.2. Reionization 


Since polarization directly probes the last scatter¬ 
ing epoch the first thing we learn is when that oc¬ 
curred, i.e. what fraction of photons last scattered at 
2 « 1000 when the universe recombined, and what frac¬ 
tion rescattered when the intergalactic medium reion¬ 
ized at z r i > 5. 


Since rescattering erases fluctuations belo w the hori- 
zon scale and regenerates them only weakly ( Efstathior] 
1988), we already know from the reported excess (over 
COBE ) of sub-degree scale anisotropy that the optical 
depth during the reionized epoch was r < 1 and hence 


2ri < 100 


n 0 h 2 \ 1/3 
0.25 J 


n B h 2 y 2/3 
0.0125 J 


( 11 ) 


The main reason why polarization is so useful to the 
reconstruction process is that (cosmologically) it can 
only be generated by Thomson scattering. The polar¬ 
ization spectrum of the CMB is thus a direct snapshot of 
conditions on the last scattering surface. Contrast this 
with temperature fluctuations which can be generated 
by changes in the metric fluctuations between last scat¬ 
tering and the present, such as those created by gravi¬ 
tational potential evolution. This is the reason why the 
mere detection of large-angle anisotropies by COBE did 
not rule out wide classes of models such as cosmological 
defects. To use the temperature anisotropies for the re¬ 
construction problem, one must isolate features in the 
spectrum which can be associated with last scattering, 
or more generally, the universe at a known redshift. 

Furthermore, the polarization spectrum has poten¬ 
tial advantages even for extracting information from 
features that are also present in the temperature spec¬ 
trum, e.g. the acoustic peaks. The polarization spec¬ 
trum is generated by local quadrupole anisotropies 
alone whereas the temperature spectrum has compa¬ 
rable contributions from the local monopole and dipole 
as well as possible contributions between last scattering 
and the present. This property enhances the promi¬ 
nence of features (see Fig. 0 ) as does the fact that 
the scalar polarization has a relatively sharp projection 
due to the geometry (see §[h^). Unfortunately, these 
considerations are mitigated by the fact that the polar¬ 
ization amplitude is so much weaker than the tempera¬ 
ture. With present day detectors, one needs to measure 
it i n broad l-bands to increase the signal-to-noise isee 

Fi eiy and Tabic |). 


It is thus likely that our universe has the interesting 
property that both the recombination and reionization 
epoch are observable in the temperature and polariza¬ 
tion spectrum. 

Unfortunately for the temperature spectrum, at these 
low optical depths the main effect of reionization is an 
erasure of the primary anisotropies (from recombina¬ 
tion) as e~ T . This occurs below the horizon at last 
scattering, since only on these scales has there been 
sufficient time to convert the originally isotropic tem¬ 
perature fluctuations into anisotropies. The uniform 
reduction of power at small scales has the same effect 
as a change in the overall normalization. For z T i ~ 5- 
20 the difference in the power spectrum is confined to 
large angles {l < 30). Here the observations are lim¬ 
ited by “cosmic variance”: the fact that we only have 
one sample of the sky and hence only 21 + 1 samples 
of any given multipole. Cosmic variance is the dom¬ 
inant source of uncertainty on the low-1 1 temperature 


spectrum in Fig. 14 


The same is not true for the polarization. As we have 
seen the polarization spectrum is very sensitive to the 
epoch of last scattering. More specifically, the location 
of its peak depends on the horizon size at last scattering 
and its height depends on the duration of last scattering 
(Efstathiou 1988). This signature is not cosmic variance 
limited until quite late reionization, though the combi¬ 
nation of low optical depth and partial polarization will 
make it difficult to measure in practice (see Fig. 0 and 
Table 0 ). Optical depths of a few percent are poten¬ 


tially observable from the MAP satellite (Zaldarriaga 


et al.| 1997) and of order unity from the POLAR exper¬ 


iment (Keating et al. 1997). 


On the other hand, these considerations imply that 
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the interesting polarization signatures from recombina¬ 
tion will not be completely obscured by reionization. 
We turn to these now. 


4.3. Scalars, Vectors, &: Tensors 


There are three types of fluctuations: scalars, vec¬ 
tors and tensors, and four observables: the temperature, 
B-mode, B-mode, and temperature cross polarization 
power spectra. The CMB thus provides sufficient infor¬ 
mation to separate these contributions, which in turn 
can tell us about the generation mechanism for fluctu¬ 
ations in the early universe (see §p~f|). 


Ignoring for the moment the question of foregrounds, 
to which we turn in §5.2, if the B-mode polarization 


greatly exceeds the B-mode then scalar fluctuations 
dominate the anisotropy. Conversely if the B-mode is 
greater than the B-mode, then vectors dominate. If 
tensors dominate, then the B and B are comparable 
(see Fig. ID- These statements are independent of the 
dynamics and underlying spectrum of the perturbations 
themselves. 


The causal constraint on the generation of a quad- 
rupole moment (and hence the polarization) introduces 
further distinctions. It tells us that the polarization 
peaks around the scale the horizon subtends at last 
scattering. This is about a degree in a flat universe 
and scales with the angular diameter distance to last 
scattering. Geometric projection tells us that the low-£ 
tails of the polarization can fall no faster than £ e , £ 4 
and & 2 for scalars, vectors and tensors (see §[c|). The 
cross spectrum falls no more rapidly than £ 4 for each. 
We shall see below that these are the predicted slopes 
of an isocurvature model. 


Furthermore, causality sets the scale that separates 
the large and small angle temperature-polarization cor¬ 
relation pattern. Well above this scale, scalar and vec¬ 
tor fluctuations should show anticorrelation (tangential 
around hot spots) whereas tensor perturbations should 
show correlations (radial around hot spots). 

Of course one must use measurements at small enough 
angular scales that reionization is not a source of confu¬ 
sion and one must understand the contamination from 
foregrounds extremely well. The latter is especially true 
at large angles where the polarization amplitude de¬ 
creases rapidly. 


4.4. Adiabatic vs. Isocurvature Perturbations 

The scalar component is interesting to isolate since 
it alone is responsible for large scale structure forma¬ 


tion. There remain however two possibilities. Density 
fluctuations could be present initially. This represents 
the adiabatic mode. Alternately, they can be generated 
from stresses in the matter which causally push matter 
around. This represents the isocurvature mode. 


The presence or absence of density perturbations 
above the horizon at last scattering is crucial for the fea¬ 
tures in both the temperature and polarization power 
spectrum. As we have seen in §3.3.1, it has as strong ef¬ 
fect on the phase of the acoustic oscillation. In a typical 
isocurvature model, the phase is delayed by tt/2 mov¬ 
ing structure in the temperature spectrum to smaller 
angles (see Fig. |l^). Consistency checks exist in the 
B-polarization and cross spectrum which should be out 
of phase with the temperature spectrum and oscillat¬ 
ing at twice the frequency of the temperature respec¬ 
tively. However, if the stresses are set up sufficiently 
carefully, this acous tic pha se test can be evaded by an 
isocurvature model (Furok 1996). Perhaps more impor¬ 
tantly, acoustic features can be washed out in isocurva¬ 
ture models with complicated small scale dynamics to 
force the acoustic oscillation, as in many defect models 
( Albrecht et al. 1996). 

Even in these cases, the polarization carries a robust 
signature of isocurvature fluctuations. The polarization 
isolates the last scattering surface and eliminates any 
source of confusion from the epoch between last scat¬ 
tering and the present. In particular, the delayed gen¬ 
eration of density perturbations in these models implies 
a steep decline in the polarization above the angle sub¬ 
tended by the horizon at last scattering. The polar¬ 
ization power thus hits the asymptotic limits given in 
the previous section in contrast to the adiabatic power 
spectra shown in Fig. 0 which have more power at large 
angles. To be more specific, if the B-power spectrum 
falls off as l 6 and or the cross spectrum as f? 4 , then the 
initial fluctuations are isocurvature in nature (or an adi¬ 
abatic model with a large spectral index, which is highly 
constrained by COBB). 

Of course measuring a steeply falling spectrum is dif¬ 
ficult in practice. Perhaps more easily measured is the 
first feature in the adiabatic B-polarization spectrum at 
twice the angular scale of the first temperature peak and 
the first sign crossing of the correlation at even larger 
scales. Since isocurvature models must be pushed to 
the causal limit to generate these features, their obser¬ 
vation would provide good evide nce for the adiabatic 
nature of the initial fluctuations ( Hu, Spergel fe White 
1997). 
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Fig. 15.— Tensor power spectra. The tensor temperature, temperature-cross-polarization, 15-mode, and H-mode 
polarization for a scale invariant intial spectrum of tensors with Ho = 1.0, £Ib = 0.1 and h. = 0.5. The normalization 
has been artificially set (high) to have the same quadrupole as the scalar spectrum of Fig. [hi]. Also shown is the 
MAP lcr upper limits on the U-rnode due to noise, i.e. assuming no signal, showing that higher sensitivity will 
likely be necessary to obtain a significant detection of B-polarization from tensors. 
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4.5. Inflation vs. Defects 

One would like to know not only the nature of the 
fluctuations, but also the means by which they are gen¬ 
erated. We assume of course that they are not merely 
placed by fiat in the initial conditions. Let us first divide 
the possibilities into broad classes. In fact, the distinc¬ 
tion between isocurvature and adiabatic fluctuations is 
operationally the same as the distinction between con¬ 
ventional causal sources (e.g. defects) and those gener¬ 
ated by a period of superluminal expansion in the early 
universe (i.e. inflation). It can be shown that inflation 


is tjhe only causal mechanism for generating supcilioii 


zon siz e d e nsity (curvature) fluctua t ions (Liddk:| 1995) 


Since the slope of the power spectrum in E can be traced 
directly to the presence of “superhorizon size” temper¬ 
ature, and hence curvature, fluctuations at last scat¬ 
tering, it represents a “test” of inflation ( Hu fe White| 
1997, Bpergel & Zaldarriaga 1997). The acoustic phase 
test, either in the temperature or polarization, repre¬ 
sents a marginally less robust test that should be easily 
observable if the former fails to be. 

These tests, while interesting, do not tell us anything 
about the detailed physics that generates the fluctu¬ 
ations. Once a distinction is made between the two 
possibilities one would like to learn about the mecha¬ 
nism for generating the fluctuations in more detail. For 
example in the inflationary case there is a well known 
test of single-field slow-roll inflation which can be im¬ 
proved by using polarization information. In principle, 
inflation generates both scalar and tensor anisotropies. 
If we assume that the two spectra come from a single 
underlying inflationary potential their amplitudes and 
slopes are not independent. This leads to an algebraic 
consistency relation between the ratio of the tensor and 
scalar perturbation spectra and the tensor spectral in¬ 
dex. However information on the tensor contribution to 
the spectrum is limited by cosmic variance and is easily 
confused with other effects such as those of a cosmolog¬ 
ical constant or tilt of the initial spectrum. By using 
polarization information much smaller ratios of tensor 
to scalar perturbations may be probed, with more ac¬ 
curacy, and the test refined (see Zaldarriaga et al. 1997 
for details). In principle, this extra information may 
also allow one to reconstruct the low order derivatives 
of the inflaton potential. 

Similar considerations apply to causal generation of 
fluctuations without inflation. Two possibilities are a 
model with “passive evolution” where initial stress fluc¬ 
tuations move matter around and “active evolution” 


where exotic but causal physics continually generates 
stress-energy perturbations inside the horizon. All cos¬ 
mological defect models bear aspects of the latter class. 
The hallmark of an active generation mechanism is the 
presence of vector modes. Vector modes decay with the 
expansion so in a passive model they would no longer 
be present by recombination. Detection of (cosmologi¬ 
cal) vector modes would be strong evidence for defect 
models. Polarization is useful since it provides a unique 
signature of vector modes in the dominance of the B- 
mode polarization ( Hu fe White| 1997). Its level for 
severa l specific defect models is given in |Seljak, Pen 


Turok| (1997). 


5. Phenomenology 
5.1. Observations 

While the theoretical case for observing polariza¬ 
tion is strong, it is a difficult experimental task to ob¬ 
serve signals of the low level of several fiK and be¬ 
low. Nonetheless, polarization experiments have one 
potential advantage over temperature anisotropy exper¬ 
iments. They can reduce atmospheric emission effects 
by differencing the polarization states on the same patch 
of sky instead of physically chopping between different 
angles on the sky since atmos pher ic emission is thought 
to be nearly unpolarized (see §[h^). However, to be suc¬ 
cessful an experiment must overcome a number of sys¬ 


tematic effects, many of which are discussed in Keating 
|et al. (1997). It must at least balance the sensitivity of 


the instrument to the orthogonal polarization channels 
(including the far side lobes) to nearly 8 orders of mag¬ 
nitude. Multiple levels of switching and a very careful 
design are minimum requirements. 

To date the experimental upper limits on polariza¬ 
tion of the CMB have been at least an order of mag¬ 
nitude larger than the theoretical expectations. The 
original polarization limits go back to Penzias & Wil¬ 
son (1965) who set a limit of 10% on the polarization 
of the CMB. There have been several subsequent up¬ 
per limits which have now reached the level of ~ 20/iK 
(see Tab. |l!) , about a factor of 5-10 above the predicted 
levels for popular models (see Tab. j^). 

5.2. Foregrounds 

Given that the amplitude of the polarization is so 
small the question of foregrounds is even more impor¬ 
tant than for the temperature anisotropy. Unfortu¬ 
nately, the level and structure of the various foreground 
polarization in the CMB frequency bands is currently 
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Reference 

Frequency (GHz) 

Ang. Scale 

AT(/xK) 

(CL) 

Penzias & Wilson (1965) 

4.0 

- 

To 5 

95% 

Nanos (1979) 

9.3 

15° beam 

2 x 10 3 

90% 

Caderni et al. (1978) 

100-600 

0.5° < 9 < 40° 

2 x 10 3 (40°/(9) 1/2 

65% 

Lubin & Smoot (1981) 

33 

£<, 3 

200 

95% 

Partridge et al, (1988) 

5 

V <9<3' 

100 

95% 

Wollack et al. (1993) 

26-36 

50 < £ < 100 

25 

95% 

Netterfield et al. (1995) 

26-46 

50 < £ < 100 

16 

95% 


Table 1: Experimental upper limits on the polarization of the CMB. The values are converted from the quoted 
A T/T by multiplying by T = 2.728K, thus these numbers only approximate the bandpowers used in Tab. |^. 


Cosmological Parameters 
n 0 n B z Ti T/S 

£: 

E 

2-20 (/xK) 

TE B 

£: 100-200 (/xK) 

E TE B 

£: 700-1000 (/xK) 

E TE B 

1.0 

0.05 

0 

0 

0.02 

(—)0.42 

0 

0.68 

(+)3.6 

0 

3.7 

(+)4.8 

0 

1.0 

0.10 

0 

0 

0.017 

(—)0.40 

0 

0.71 

(+)4.2 

0 

3.4 

(+)3.3 

0 

1.0 

0.05 

5 

0 

0.03 

(—)0.59 

0 

0.68 

(+)3.6 

0 

3.7 

(+)4.8 

0 

1.0 

0.05 

10 

0 

0.06 

(—)0.78 

0 

0.68 

(+)3.6 

0 

3.6 

(+)4-8 

0 

1.0 

0.05 

20 

0 

0.14 

(-)l.l 

0 

0.66 

(+)3-5 

0 

3.5 

(+)4.7 

0 

1.0 

0.05 

0 

0.1 

0.02 

(—)0.40 

0.01 

0.65 

(+)3-4 

0.05 

3.5 

(+)4.6 

0.01 

0.4a: 

0.05 

0.1 

0 

0.01 

(—)0.18 

0 

0.44 

(+)2-3 

0 

2.8 

(+)5-5 

0 

0.4 a 

0.05 

0 

0 

0.02 

(—)0.37 

0 

0.65 

(+)4.0 

0 

4.1 

(+)6-4 

0 


Table 2: Bandpowers for a selection of cosmological models all with h = 0.5 and scale invariant spectra save for 
the T/S = 0.1 model where T/S = 7(1 — n s ). The low-O 0 models are indicated as 0.4/< for the open model and 
0.4a for the model with = 1 — Oq = 0.6. The bandpower is defined in units of the COBE quadrupole over the 
range of £ indicated. The sign of the cross correlation is indicated in parentheses and note that the band averaged 
signal may be suppressed from cancellation of correlated and anticorrelated angular regimes. 


not well known. We review some of the observations 
in the adjacent radio and IR bands (a more complete 


discussion can be found in [Keating et alj 1997). At- 


ized, however Hildebrand fe Dragovan (1995) find that 
in their observations the majority of dust is polarized 
at the « 2% level at 100/xm with a small fraction of 


mospheric emission is believed to be negligibly polar- regions approaching 10% polarization. Moreover Keat 


ize< l_( Keating et al. 1997), leaving the main astrophysi- 
cal foregrounds: free-free, synchrotron, dust, and point 
source emissions. Of these the most important fore¬ 
ground is synchrotron emission. 

Free-free emission (bremsstrah lung) is intrinsically 
unpolarized ( Rybicki fe Lightman 1979) but can be par¬ 
tially polarized by Thomson scattering within the HII 
region. This small effect is not expected to polarize the 
emission by more than 10% (Keating et al. 1997). The 
emission is larger at low frequencies but is not expected 
to dominate the polarization at any frequency. 

The polarization of dust is not well known. In princi¬ 
ple, emission from dust particles could be highly polar¬ 


ing et al. (1997) show that even at 100% polarization, 


extrapolation of the IRAS 100/xm map with the COBE 
FIRAS index shows that dust emission is negligible be¬ 
low 80GHz. At higher frequencies it will become the 
dominant foreground. 

Radio point sources are polarized due to synchrotron 
emission at < 20% level. For large angle experiments, 
the random contribution from point sources will con¬ 
tribute negligibly, but may be of more concern for the 
upcoming satellite missions. 

Galactic synchrotron emission is the major concern. 
It is potentially highly polarized with the fraction de¬ 
pendent on the spectral index and depolarization from 
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Faraday rotation and non-uniform magnetic fields. The 
level of polarization is expected to lie between 10%-75% 
of a total intensity which itself is approximately 50/iK 
at 30GHz. This estimate follows from extrapolating the 
Brouw fe Spoelstra| (1976) measurements at 1411 MHz 
with an index of T oc v~ 3 . 

Due to their different spectral indices, the minimum 
in the foreground polarization, like the temperature, lies 
near 100GHz. For full sky measurements, since syn¬ 
chrotron emission is more highly polarized than dust, 
the optimum frequency at which to measure intrin¬ 
sic (CMB) polarization is slightly higher than for the 
anisotropy. Over small regions of the sky where one or 
the other of the foregrounds is known a priori to be 
absent the optimum frequency would clearly be differ¬ 
ent. However as with anisotropy measurements, with 
multifrequency coverage, polarized foregrounds can be 
removed. 


It is also interesting to consider whether the spatial 
as well as frequency signature of the polarization can 
be used to separate foregrounds. Using angular power 
spectra for the spatial properties of the foregrounds 
is a simple generalization of methods already used in 
anisotropy work. For instance, in the case of synchotron 
emission, if the spatial correlation in the polarization 
follows that of the temperature itself, the relative con¬ 
tamination will decrease on smaller angular scales due 
to its diffuse nature. Furthermore the peak of the cos¬ 
mic signal in polarization occurs at even smaller angular 
scales than for the anisotropy. 

One could attempt to exploit the additional prop¬ 
erties of polarization, such as its E- and E-mode na¬ 
ture. We mentioned above that in a wide class of 
models where scalars dominate on small angular scales, 
the polarization is predicted to be dominantly .E-mode 
(Kamionkowski et al. 1997 Zaldarriaga fe Seljak| 1997). 
seljak (1997) suggests that one could use this to help 
eliminate foreground contamination by “vetoing” on ar¬ 
eas of E-mode signal. However in general one does not 
expect that the foregrounds will have equal E- and E- 
mode contribution, so while this extra information is 
valuable, its use as a foreground monitor can be com¬ 
promised in certain circumstances. Specifically, if a cor¬ 
relation exists between the direction of polarization and 
the rate of change (curvature) of its amplitude, the fore¬ 
ground will populate the two modes unequally. Two 
simple examples: either radial or tangential polariza¬ 
tion around a source with the amplitude of the polar¬ 
ization dropping off with radius, or polarization parallel 
or perpendicular to a “jet” whose amplitude drops along 


the jet axis. Both examples would give predominantly 
E-mode polarization. 


5.3. Data Analysis 


Several authors have addressed the question of the 
optimal estimators of the polarization power spectra 
from high sensitivity, all-sky maps of the polarization. 
They suggest that one calculate the coefficients of the 
expansion in spin-2 spherical harmonics and then form 
quadratic estimators of the power spectrum, as the av¬ 
erage of the squares of the coefficients over m, corrected 


for noise bias as in the example of Fig. 14 


We shall return to consider this below, however be¬ 
fore such maps are obtained we would like to know how 
to analyze ground based polarization data. These data 
are likely to consist of Q and U measurements from 
tens or perhaps hundreds of pointings, convolved with 
an approximately gaussian beam on some angular scale. 
How can we use this data to provide constraints or mea¬ 
surements of the electric and magnetic power spectra, 
presumably averaged across bands in £? 

For a small number of points, the simplest and most 
powerful way to obtain the power spectrum is to per¬ 
form a likelihood analysis of the data. The likelihood 
function encodes all of the information in the mea¬ 
surement and can be modified to correctly account for 
non-uniform noise, sky coverage, foreground subtrac¬ 
tion and correlations between measurements. Opera¬ 
tionally, one computes the probability of obtaining the 
measured points Qi = Q(n,) and Ui = U(hi) assum¬ 
ing a given “theory” (including a model for foregrounds 
and detector noise) and maximizes the likelihood over 
the theories. For our purposes, the theories could be 
given simply by the polarization bandpowers in E and 
B for example, or could be a more “realistic” model 
such as CDM with a given reionization history. The 
confidence levels on the parameters are obtained as mo¬ 
ments of the likelihood function in the usual way. Such 
an approach also allows one to generalize the analysis 
to include temperature information (for the cross cor¬ 
relation) if it becomes available. 

Assuming that the fluctuations are gaussian, the like¬ 
lihood function is given in terms of the data and the 
correlation function of Q and U for any pair of the n 
data points. The calculation of this correlation function 


is straightforward, and Kamionkowski et al. (1997) dis¬ 
cuss the problem extensively. Let us assume that we 
are fitting only one component or have only one fre¬ 
quency channel. The generalization to multiple frequen- 
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cies with a model for the foreground is also straightfor¬ 
ward. We shall also assume for notational simplicity 
that we are fitting only to polarization data, though 
again the generalization to include temperature data is 
straightforward. The construction is as follows. We 
define a data vector which contains the Q and U infor¬ 
mation referenced to a particular coordinate system (in 
principle this coordinate system could change between 
different subsets of the data). Call this data vector 


from satellite experiments. For these experiments, the 
large volume of data is however an issue in the analysis 
pipeline design, as has been addressed by several au¬ 
thors. The generalization of the above analysis proce¬ 
dure to include filtering and compression is straightfor¬ 
ward, and directly analogous to the case of anisotropy, 
so we will not discuss it explicitly here. 

6. Future Prospects 


D = (Q 1 ,U 1 ,...,Q n ,U n ), (12) 


which has N = 2n components. We can construct the 
likelihood of obtaining the data given a theory once we 
know the correlation matrix Cij: 


£(D\T) oc 


1 


Vdet C 


exp 


-DjCrjDj 


(13) 


All of the theory information is encoded in Cij = Cf} + 
Nij , where Njj is the noise correlation matrix, to be 
provided by the experiment, and Cf} is a function of 
the theory parameters. 

All that remains is to compute each element of Cf } 
for a given theory. Consider a pair of points i and j 
corresponding to 4 entries of our data vector Dj. Fol¬ 
lowing Kamionkowski et al. (1997), define Q' and U' 
as the components of the polarization in a new coor¬ 
dinate system, where the great arc connecting ni and 
n 2 runs along the equator. Expressions for {Q'Q') and 
{U'U') in terms of the E and B angular power follow 
directly from the definitions of these spectra and the 
spin-weighted spherical harmonics. They can be found 


in Kamionkowski et al. (1997) [their Eqs. (5.9,5.10)]. 
They also give the flat sky limit of these equations. 
Knowing the angle faj about n, through which we must 
rotate our primed coordinate system to return to the 
system in which our data is defined we can write 


Qi = Q[ cos 2 <f>ij + U[ sin 2^ 

Ui = sin 2cj>ij — U[ cos 2(j>ij (14) 


and similarly for the jth element. Thus using the 
known expressions for {Q'Q') and {U'U') we can cal¬ 
culate ( QiQj ), ( UiUj ), {QiUj) and ( UiQj ) for the zth 
and jth pixel and thus all of the elements of Cfy. Sub¬ 
stitution of into Eq. ( |l3| ) allows one to obtain limits 
on any theory given the data. 

Finally let us note that the method outlined here is 
completely general and thus can be applied also the 
the high-sensitivity, all-sky maps which would result 


Following a hiatus of some years, experimental ef¬ 
forts to detect CMB polarization are now underway. 
At the l arge angular sc ale, an experiment based in Wis¬ 
consin ( [Keating et al. 1997) plans to make a 30GHz 
measurement of polarization on 7° angular scales with 
a sensitivity of a few /iK per pixel. The main goal of 
the experiment would be to look for the signature of 
reionization (see Fig. |T~i| j . Several groups are currently 
considering measuring the signature at smaller scales 
from recombination. 


The MAP and Planck missions also plan to mea¬ 
sure polarization. Because of their all sky nature these 
missions will be the first to be able to measure the po¬ 
larization power spectrum and temperature cross cor¬ 
relation with reasonable spectral sensitivity. We show 
an example of the sensitivity to the which MAP should 
nominally obtain in the absence of foregrounds and sys¬ 
tematic effects in Fig. |l4|. As can be seen, MAP should 
certainly detect the E-mode of polarization at medium 
angular scales, and obtain a significant detection of 
the temperature polarization cross-correlation. These 
signatures can be useful for separating adiabatic and 
isocurvature models of structure formation as we have 
seen in §|[ 

However the detailed study of the polarization power 
spectrum will require the improved sensitivity and ex¬ 
panded frequency coverage of Planck for detailed fea¬ 
tures and foreground removal respectively. Planck will 
have the sensitivity to make a measurement of the large- 
angle polarization predicted in CDM models, regardless 
of the epoch of reionization if foreground contamination 
can be removed. Likewise, it can potentially measure 
small levels of E-mode polarization. The separation of 
the E and B modes is of course crucial for the isolation 
of scalar, vector and tensor modes and so the recon¬ 
struction problem in general. 

Clearly, the polarization spectrum represents another 
gold mine of information in the CMB. Though signifi¬ 
cant challenges will have to be overcome, the prospects 
for its detection are bright. 
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